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Abstract 

In this work we present a set of microscopic U{\)y. U{\) models which realize insulating phases 
with a quantized Hall conductivity a^y The models are defined in terms of physical degrees of 
freedom, and can be realized by local Hamiltonians. For one set of these models, we find that 
(Jxy is quantized to be an even integer The origin of this effect is a condensation of objects made 
up of bosons of one species bound to a single vortex of the other species. For other models, 
the Hall conductivity can be quantized as a rational number times two. For these systems, the 
condensed objects contain bosons of one species bound to multiple vortices of the other species. 
These systems have excitations carrying fractional charges and non-trivial mutual statistics. We 
present sign-free reformulations of these models which can be studied in Monte Carlo, and we 
use such reformulations to numerically detect a gapless boundary between the quantum Hall and 
trivial insulator states. We also present the broader phase diagrams of the models. 
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1. Introduction 

A major project in condensed matter research is the discovery and classification of new topo- 



■oje 

logical phases. ill 4l3ll One class of these phases are the topological insulators. Though much 
progress has been made in the study of topological insulators of free fermions,|[ll|2| less is known 
about interacting systems. Much work has focused on understanding topological insulator phases 
that do not have long-range entanglement, also known as Symmetry Protected Topological (SPT) 
phases, i^l-l^ and very recent studies have also pursued cases with long-range entanglement. lfl4l- 
1611 . Recent works 191 1 17141911 have used Chern-Simons approaches to provide understanding of 
interacting topological insulators in {2+VjD, and in particular of so-called integer Quantum Hall 
states of bosons (for a review, see Ref. ll20ll ). Several papers have proposed qualitative construc- 



tion of such phases using Chern-Simons flux attachment 1 1711 and slave-particle approaches. 1121 



2311 However, so far there have been no microscopic models producing such states. 

Here we will present fully tractable realizations of both integer and fractional Quantum Hall 
phases with U{1) x U{1) symmetry. Our approach is different from Chern-Simons and slave- 
particle approaches, in that we work directly with physical degrees of freedom and do not intro- 
duce artificial fluxes or enlarge the Hilbert space. More specifically, we think directly in terms of 
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charge and vortex degrees of freedom, which are all precisely mathematically defined in our lat- 
tice models. In this way, our work is close in spirit to pursuits to understand fractionalized phases 



of spins and bosons in terms of the vortex phvsics. rt24i 12511 The U{1) x U{1) structure allows 
us to provide an unambiguous and simple physical picture of the integer and fractional quantum 
Hall states of bosons. Thus, an elementary integer quantum Hall state is obtained when a vortex 
in one species binds a charge of the other species and the resulting composite object condenses; 
our approach provides a precise meaning of such a condensation. General integer quantum Hall 
states are obtained when a vortex in one species binds a fixed number of charges of the other 
species. On the other hand, fractional quantum Hall states are obtained when we have conden- 
sation of composite objects that are bound states of d vortices of one species and c particles of 
the other species. In this case we show that the system has a fractional quantum Hall response 
given by a^y = 2c/ d and has quasiparticles carrying fractional charges of l/d of the micro- 
scopic charges and non-trivial mutual statistics. Such d-tupled vortex condensation leading to 



charge fractionalization is reminiscent of the idea in Refs. II241I25I1 of paired- vortex condensation 
leading to Z2 fractionalized phases, although in the present case there is also binding of charges 
and the phase shows quantized Hall response. The integer quantum Hall phases discussed in 
this work are examples of SPT phases, while the fractional quantum Hall states are examples of 
interacting topological insulators with topological order and long-range entanglement. 

An interesting aspect of our models is that they are examples of interacting topological in- 
sulators that can be studied in a sign-free Monte Carlo, thus adding to a growing list of exotic 



quantum phases and transitions that can be explored numeric allv. 112611 Furthermore, perhaps for 
the first time, we can simulate systems that have non-trivial quantum Hall responses and in partic- 
ular gapless edge states. In this paper, we perform such a study of a boundary between Quantum 
Hall and trivial insulators and provide compelling evidence of gaplessness at the boundary. 

The paper is organized as follows. In Sec. |2]we describe our models and demonstrate that 
they realize integer and fractional quantum Hall phases. In Sec. |3]we derive the sign-free refor- 
mulations used in the numerics. In Sec. |4]we present the results of these numerics, including 
evidence for gapless modes at the boundary of the quantum Hall states. In Sec.|5]we show the 



broader phase diagrams for our models. Appendix A [reviews a duality transformation which is 



used throughout this paper In [Appendix B we demonstrate the connections between our mod- 
els and various field-theoretic approaches, specifically non-linear sigma models with topological 



terms and A'-matrix theories. Finally, in Appendix C we present a local Hamiltonian which 
gives the action discussed in the main text. 

2. Explicit Models with Integer and Fractional Quantum Hall Effect 

In this section we introduce the models which are the subject of this work. We then rewrite the 
models in terms of new variables, in the process demonstrating that the models realize quantum 
Hall phases. The models are given by the following action in (2h-1)D Euclidean space-time: 



r,r' R,R' 

+ iY,w{R-R')[V -Kji]{R)-J2{R')- (1) 



R.R' 



The Ji variables are integer-valued conserved currents residing on a cubic lattice (the 'direct' 
lattice) labeled with the index r. They satisfy V • J7i =0, which means that they form closed 
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loops. The J2 variables are the same except that they reside on a lattice dual to the direct lattice, 
and whose sites are labeled by the index R. These variables represent space-time currents of two 
separately conserved species of bosons. In (2+l)D, applying the curl operator to a vector object 
on the direct lattice, such as a J\ (r) variable, leads to a vector object on the dual lattice, denoted 
[Vx 

The first two terms in the above action are short-ranged repulsive interactions within each 
species of boson. The third term is a short-ranged interaction between currents of different 
species. The latter term may appear unusual, but this action is local, and in Appendix C we 



show that it can arise from a local Hamiltonian. Also, even though it appears that this term 
treats the J\ and J2 currents differently, we can integrate by parts and see that, as long as 
w(R, R') = w{R - R') = w{R' - R), the Ji and J2 enter symmetrically. 

In loop actions such as Eq. ([T]i, it is possible to make a change of variables to obtain equivalent 
actions that may be easier to interpret. One such change of variables is a well-known duality 



transformation, (given in Appendix A 1 which can be generalized l34ll to arbitrary modular group 
transformations. 1 35I 36 1 This is the change of variables we will use to show that our model can 
produce integer and fractional bosonic quantum Hall insulators. 

To make this change of variables, we first rewrite Eq. dl] in fc-space as follows: 



S = 



hWl^i(^')l'+^'2WI^2(fc)|'] +iY.e{k)M-k)-dj2{k). (2) 



In the above expression, Ji(fc) = ■Jiij)e.~^^'^ / y/Yol and wi(fc) = ^^vi{r—r')e~^^''^'^~'^'\ 
J2{k) and V2{k) are defined similarly to Ji{k) and wi(fc), but on the dual lattice indexed by R. 
The "gauge" field dj2{r) is defined such that J2{R) = [V x dj2]{R) (such representation 
is possible given the divergenceless conditions and also conditions of the vanishing total cur- 
rents discussed below, and the action is independent of the gauge choice ll33[|37ll ). Starting from 
Eq. O, we obtain e{k) = \ fk?Y.Rw{R - i?')e~*'^"^'\ where l/^p = ^^(2 - 2cosfc^) 
arises from the lattice derivatives and the sum is over all directions. We assume that vi (fc), V2 (fc) 
and 6{k) are even in fc. 

In order to determine the physical properties of the model, we couple it to external gauge 
fields by adding the following terms to Eq. 



6S 



AT\k) + J2{~k) 



^cxt 



(fc-) 



(3) 



We can use the duality transform from Appendix A to go from the Ji variables to dual Qi 
variables as follows: 



2^Qi(fc) + 9{k)J2{k) + [V X If t](fc) 



2 ^ 

k 



\fk\'v,{k) 



2^ 

fc 



(fc) IJ2 (fc) I ' + * E ^2 (" fc) • (fc) ■ 



(4) 



If we interpret the Ji variables as physical bosons, the Qi variables can be interpreted as vortices 
in the boson phase variables. Like the J'l variables, the Qi variables are divergenceless and 
therefore form closed loops. A technical remark: We have found it convenient to require zero 
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total current in our system, Jitot = J2r ~ ^ and J2tot = 0. This makes the above duality 

procedure exact, |33l 37 1 and requires Qitot = 0. All of the currents defined below will also 
satisfy this condition. 



As discussed in our previous work. 11341 in the above action we can make the following change 
of variables: 



02 = cQi~dj2. 

This change of variables is valid if the matrix 

e PSL{2,Z), 



(5) 
(6) 



(7) 



i.e., a, b, c, d are integers such that ad — be ^ 1. Since the above matrix is an element of the 
modular group, we call this change of variables a modular transformation and will often refer to it 
simply (a, 6, c, d). We can then perform the duality transform to go from the Fi variables to dual 
Qi variables, which gives us an action in terms of the Qi and Q2 variables. This transformation, 
from Ji, J2 variables to Qi, Q2 variables, is the generalization of the duality operation to modular 
transformations. After performing this change we are left with the following action: 



S = 



k 



k 



Vgi{k) 



Gi{k) 



c[V X Af*](fc) 



27r 



Q2{k) 



c[V X If *](fc) 



27r 



'e(fc)^l(-fc)•ae2(fc)-^^ 

k 

eg{k)c 



■:[2Tra- 0g{k)c 



■[V X •If'(fc) 



a 



27r 



where C/2 = V x dg2 and 

vgi/2{k) 

27r 



g^{-k) ■ AT\k) + g2{-k) ■ A-'^^k) 



{2Tl)\,^,2{k) 



[2TTd + e{k)cY + 

[2TTb + e{k)a] [2nd- 



Vl{k)v2(k)\fk\^C^ 

h9{k)c]+vi{k)v2{k)\fk\''ca 



(8) 



(9) 



[2Trd + e{k)c]^ + vi{k)v2{k)\fk\^c^ 



We can deduce the quantum Hall properties of the system from the action in terms of the 
Q variables. If we ignore the external fields Ay^^, we see that this action has the same form as 
Eq. (|2]i, but the action is defined in terms of different variables which see different inter-particle 
potentials. If vg is large, so that only small loops in the G variables can form, we say that 
the Q variables are 'gapped', and in this phase we can interpret the Q variables as describing 
gapped quasiparticles above the ground state. Equation (|8) is the action for these quasiparticles. 
Restoring the external fields j4^y2 ^"d integrating out the gapped Q variables, we would only 
generate Maxwell-like terms for the external fields (assuming that vg potentials are short-ranged, 
which is true for d ^ 0). Hence, we can extract Hall-like conductivity describing transverse-cross 
response from the fourth term in the above action: 



rllik) 



2c[2Tra- 0g{k)c] 
2^ 



(10) 



This conductivity is defined in units of /h. 

Now consider the conductivity in the A; — > limit. Here and below we assume d 7^ 0. (If we 
use the above transformation with d = 0, the Q particles are essentially vortices in the J, and if 
they are gapped, this gives a superfluid and will be discussed in Sec.|5]) From Eq. (|9]l we can see 
that in this limit and for short-ranged v{k), 

Bg =2TT- + 0{k^). (11) 

a 

This describes the mutual statistics of the Qi and Q2 particles. We can use 9g to determine the 
conductivity in the k ^ Q limit: 

= 2^. (12) 

We see that for short-ranged potentials, we have a universal, rational conductivity. From the last 
term in Eq. (|8]l we can extract the charges of the Qi and Q2 quasiparticles to be 



charge = i (13) 
a 

relative to A^^* and Aj^* respectively. We see that when d = 1, the system has a Hall conductiv- 
ity quantized to be an even integer and excitations carrying integer charges. We propose that this 
is a realization of the bosonic integer quantum Hall effect. |@] When d > 1, we see that the Hall 
conductivity is quantized as a rational number and the Q quasiparticles carry fractional charge. 
Therefore this phase is a fractional quantum Hall effect for bosons. 

The above results are true only in the region where the Q variables are gapped. For a general 
set of potentials Vi/2{k),9{k) and coefficients a,b,c,d, we must use numerics to determine 
whether such a gapped region exists, and where it can be found. We now present a simple, short- 
ranged model in which a gapped region always exists and we know where the gapped region is 
located. This model also has the feature that the above statistical angle 9g{k) and charges of the 
quasiparticles are independent of k. This will simplify our analysis, though we expect the results 
we derive for these specific choices will hold also for any other short-ranged potentials in the 



/c — )- limit, as long as the Q particles are gapped. (In Appendix B we show how the above 
analysis can be viewed as a derivation of an effective /\ -matrix-like theory at long wavelengths 
for our quantum Hall phases.) 



To get these simple models we choose: 



^2/1 



vi/2ik) = ^^V-, (14) 

where Ai and A2 are real-valued parameters. The potentials also contain parameters c and d 
which are integer-valued. If we choose the coefficients c and d in the transformation Eq. (|6) 
[leading to Eq. ©I to be the same as those in these potentials, the resulting action for the Q 
variables will have the simplifying features described above. As we will see in a moment, in this 
case the parameters c and d give the conductivities and fractionalized charges of the quantum 
Hall phases occurring at small Ai 2, and therefore they are useful labels of the resulting quantum 
Hall phases. We will explain how these potentials were determined in Sec.|5j 
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With the above choices, the potentials in the action in terms of the Q variables take the 
following simple form: 



vg,{k) = (16) 

Og = 27r^. (17) 

a 

Therefore we can easily see that with the above choices we will always have a phase where the Qi 
variables are gapped in the limit of small A;. In this limit the system will be a bosonic quantum 
Hall insulator 

If we insert the above expressions for v{k) and d{k) into Eq. (|4]i, we see that we can write 
that action in the following way: 

k R ^ 

Here ri{R) = c/d everywhere in the system, though later we will consider spatially varying 
77(i?). Studying this action at small Ai and A2 can tell us about the physics of the quantum Hall 
states. 

First consider the situation where c = 0, which leads to i]{R) = 0, and consider the limit 
of small Ai, A2. We can see from Eq. (fT2] i that ct^^ vanishes and therefore this system is not 
in a Quantum Hall state. There is a small energy cost for loops in the Qi variables, but a large 
energy cost for loops in the J2 variables. Therefore the J2 loops are gapped. The Qi variables 
are 'condensed', which means that large loops of these variables can form. In general, since J'l 
and Qi are related by a duality transformation, if one of them is gapped the other is condensed. 
Therefore the Ji variables are gapped in this phase. Since the physical Ji and J^2 variables are 
gapped, the system is a trivial insulator This situation is shown on the right side of Fig. [T] We 
could have arrived at this conclusion also more quickly by examining Eqs. (fl4l i-(fT5Tl for c = 0, 
but it was convenient to develop the picture in the Qi and J2 variables. 

Now consider the case where c 7^ 0, so that i]{R) = c/d. Since the G variables are gapped we 
have a quantized Hall conductivity, and the system is in the bosonic quantum Hall state. Notice 
that composite objects with J2 = c and Q,\ = d see a very small energy cost, and therefore 
large loops of such objects can form. On the other hand, both the J2 and Qi variables see 
large potentials if they exist independently, so only small loops of these variables can form by 
themselves. This is illustrated on the left side of Fig. [1] 

We would like to know what carries the charges that leads to the cr^^ 7^ 0. By analogy with 
the fermionic Quantum Hall effect, we expect that the charges are being carried by edge states. 
We will examine the physics of the formation of edge states by including a boundary between 
the quantum Hall state and the trivial insulator in our system. This is accomplished by allowing 
7](i?) to vary in space. Therefore we have one region where r\(K) = c/d, and we have another 
region where 'i]{R) = 0. Note that we have defined ri{R) as varying on the dual lattice denoted 
by the index R. Before we allowed rj to vary in space, we had a symmetry between the Ji and J2 
variables in the case where vi = V2- However the different variables see the boundary differently 
and it breaks this symmetry. 

Now consider what happens at the boundary of the quantum Hall state. For example, consider 
the case where ri{R) = 1 in the region where it is non-zero. In the quantum Hall region we will 
have large loops with {J2 = l,Qi = 1), while in the trivial insulating region we have large loops 
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Figure 1: A sketch of the situation where r]{R) is allowed to vary in space, with r?{-R) = 1 in the quantum Hall region 
and 'ri{R) = in the trivial insulator region. In the trivial insulator region, loops of the Qi variables can be large, while 
the J72 variables can only foiin small loops. In the quantum Hall region, large loops of composite objects with (Qi = 1, 
J'2 = 1) exist, while loops of only Qi or only variables are small. In order to have closed loops of the J72 variables, 
there must be behavior on the edge different from that in the bulk, which leads to the gapless edge supporting non-zero 
Hall conductivity. 



of only Qi variables. This situation is illustrated in Fig.[T] The large loops of the Qi variables can 
pass through the boundary. However, the J72 loops must be bound to the Qi loops in one region 
and must disappear in the other region. The system must find a way to accomplish this while 
satisfying the constraint that all currents must be divergenceless. To do this, it exhibits behavior 
which would be energetically forbidden in the bulk. For instance, the J2 currents could run 
along the edge, as seen in Fig.[I] In this case on the right edge, the J2 currents run from places 
where the Qi variables cross the boundary from left to right to places where they cross from 
right to left. Alternatively, loops of the Q variables could be forbidden from crossing between 
the two regions, which means that vortices in the boson phase variables would be forbidden on 
the boundary. This unusual behavior leads to gapless modes on the boundary between a quantum 
Hall region and a trivial insulator In Sec. |4]we will develop mathematical description of this 
behavior and will numerically find evidence for these gapless modes. 



3. Monte Carlo Study of a Model with a Boundary 

In order to further characterize these bosonic quantum Hall phases, we now seek to find gap- 
less edge modes numerically using Monte Carlo simulations.! 38] We include edges by allowing 
r]{R) to vary in space, in one direction, which we will label as the x-direction. In our large 
L X L X L system, we will say that 
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where X is the x-coordinate on the dual lattice. The quantum Hall behavior occurs in the region 
with non-zero 77. We will take half of the system to be in the quantum Hall phase and the other 
half in the trivial insulator phase. 

Note that the original action in terms of jTi , currents, Eq. has a sign problem. On the 
other hand, the action in terms of Qi, J'2, Eq. ( fTSl l, is a statistical mechanics model that can be 
directly studied in Monte Carlo. However, we can see that the Qi variables have long-ranged 
interactions, and in systems with such interactions the time required for the simulations scales 
as L^. We would like to have only short-ranged interactions in our system, so that the Monte 
Carlo can run in a time proportional to L^, and larger sizes can be studied. The remainder of this 
section will be devoted to developing several exact reformulations of Eq. (fTSb which have only 
short-ranged interactions. These reformulations are also useful because they provide access to 
boson Green's functions. 

Consider the following action: 

S[ai,a2,Ji,J2] = i^Aiai(r)2 + i^A2a2(i?)' (20) 

r R 



27r ' 

R r R 



+ * ^ ^ [V X «i] (i?) • 52 (i?) + i ^ Ji (r) • ai (r) + * ^ ^2 (i?) • «2 (i?) • 



The new ai and a2 variables are real-valued vector fields residing on links of the direct and 
dual lattices respectively. One can check that the partition sum which results from this action 
is the same as Eq. ([T]i if we integrate out both of these species of vector fields. Note that the 
ai and a2 variables are not some gauge fields; rather, they are some local massive fields and 
are integrated over with no constraints. Appendix C| shows that such an action can be realized 



as a path integral for a local Hamiltonian with bosons coupled to oscillator degrees of freedom. 
Appendix B motivates the connection with Eq. ([T]i by starting from an action in terms of vortex 



variables, while we can also derive the result Eq. ([18) for general 'ri{R) from a direct analysis 
below. 

In order to get a sign-free action that we can study efficiently in Monte Carlo, we start with 
Eq. (ISTT i but perform different integrations. First, we integrate out the a2 variables. Then, we 
would like to integrate out the J'l variables. However, the J'l variables are integer-valued and 
constrained to be divergenceless with no total current. We enforce these constraints by adding 
new variables to our partition sum as follows. To enforce the divergenceless of the J'l variables 
we add the following term: 

/•TT 

(21) 



(We are ignoring overall constants here and below.) This introduces a 27r-periodic 0i (r) variable 
on every site of the lattice. These variables correspond to the phases of the type-1 bosons. We 
enforce the constraint that there must be no total current (in our full system) by adding another 
term to the partition sum: 

3 />7r 

'^Jit„t=o = n / djif,exp[-i-fif,'^Sr^=oJif,{r)]. (22) 

This term introduces a 27r-periodic 71^ variable for each direction /i on the lattice. This variable 
means that instead of periodic boundary conditions we have a fluctuating boundary condition 
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such that across the boundary the 4> variables differ by 7; here we chose the boundary plane 
perpendicular to the x direction to be at a; = 0, and similarly for the other directions. Note that 
these boundary terms are not related to the edges between quantum Hall and trivial insulator 
regions in the system. They are only included to make the duality procedure between the J and 
Q variables precise, so that we are confident that our various reformulations simulate exactly the 
same system for any finite size. 

Now that the Ji variables are unconstrained, we can go from integer- valued Ji to real-valued 
ji by using the following relation: 

/+00 
dji^(r) exphz2^pi^(r)ji^(r)] [...]. (23) 

This term introduces integer- valued pi^ (r) variables on every link on the lattice, and these vari- 



ables are free of any constraint. In formal duality maps. 1127143311 the physical meaning of the pi 
variables is that their curl gives vorticity in the phase variables conjugate to Ji, i.e., Qi — Vxpi. 
With all the constraints implemented, we can now integrate out the ji variables to get: 

^ = E'E ^-^i r n dj,^e-''^^^'^^'P^'^-\ (24) 

r 

+ ^ E - vim^ X p,){R) 

^ R 



(25) 



Here the primed sum over J2 variables is subject to the constraints mentioned earlier; we have 
also defined 

wi^(r) = 01 (r + A) - (f)i{r) - 7ip(5r„=o, (26) 

and used V x cDi = 0. This action is sign-free and all interactions are short-ranged, so it can be 
studied efficiently in Monte Carlo; it will allow us to detect gaplessness of the quantum Hall edge 
by looking at the spatial correlations of the 01 variables. The action is not explicitly 27r-periodic 
in the oji variables, but all physical measurements have the required periodicity. Specifically, the 
periodicity (f>i{r) — >■ (f>i{r) + 27r iV(r) and 71^ — > 71^ + 27rA/^, with integers N{r) and A'/^, can 
be accounted for by shifting the summation variables pifj.{r) Pifiir) + V^-/V(r) — Mf^Sr^^=Q, 
which does not change V x pi in the second term of Eq. ( |25] |. We will present Monte Carlo 
measurements of the e"^^ correlators at the edge in the next section. We can also use Eq. ( l25l l 
to directly show equivalence with the action Eq. (fTSl l. Indeed, if we separate the 01, pi system 
into spin-wave part and vortex part Qi = V x pi, and integrate out the spin-wave part (see e.g., 
[Appendix A| i, we would obtain exactly Eq. ( fTSl l. 

We would also like to measure Greens functions of the type-2 bosons. To get a reformulation 
which provides access to the phase variables of the type-2 bosons, (j>2, we go back to Eq. ( 1211 1 
and integrate by parts on the second line as follows: 

^77(^)[V X di](i?) • a^iR) = ^(V X [,7(i?)«2(i?)])(r) • di(r) , (27) 
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which already suggests that the two boson species will see the edge differently. We then integrate 
out the ai and J'2 variables using similar methods to those above. The resulting action is: 



—T 

2Ai ^ 



Mr) - (V X [r^{R)p2{Rmr) + ^(V x [r;(i?)^2(i?)]) W 



(28) 



where ijj2^{R) is defined similarly to Eq. ( |26] |. For the r/(i?) in Eq. ( fT9] l we can calculate: 



(Vxby(i?)^2(i?)])(r) = 3(<5^ 



(29) 



which is non-zero only near the boundaries. In the last line uj2z{R) is at the appropriate edge 
{X = Xfi or: X ~ Xl) and is calculated on the dual lattice link perpendicular to the direct lattice 
link (?', r + y) next to it, and similarly for Ld2y{R)- Thus we have the extra boundary cl>2 term in 
Eq. ( |28] |, which comes from differentiating 

Note that the 27r-periodicity in the u>2 variables is accounted for by the p2 variables. A more 
subtle observation is that in the limit of A2 = 0, the (j>2 variables have periodicity of 27r/c. 
Indeed, if we shift 02 (^) by 2ttN{R)/c with integer-valued N{R), then 0)2 (-R) is shifted by 
2ttW N (R) / c. We can simultaneously shift the summation variables P2{R) by —bWN{R) and 
observe that [a;2(^) - 27rp2(-R)]??(^)/(27r) is shifted by V N{R){ad/c)ri{R), where we used 
ad — be = 1. The latter shift is an integer-valued vector field [since {ad/c)ri{R) is an integer 
everywhere], and its curl can be absorbed into the redefinition of Ji{r), thus keeping the Ai term 
in Eq. ( |28] | unchanged. We will need to keep this in mind when measuring correlation functions 
of the (f)2 variables. 

In the case where d = 1 so that rj{R) is an integer everywhere, we can obtain one additional 
reformulation from Eq. (l2Ft . Consider making the following change of variables: 



M(r) = Mr) - (V X [T^{R)p2{R)]){r) . 



(30) 



This is an allowed change of variables since rj{R)p2{R) is an integer vector field and therefore 
its curl is a divergenceless integer field. Note that if rj{R) is not an integer, such a change of 
variables is not allowed since it will lead to a non-integer M. After making this change of 
variables, we can perform a summation over p2 to arrive at the following action: 



5[02,72,A1] = ^Fvmain[c^2M(i?); A2] + 



+00 



exp[-Vvinain(6'; A)] = ^ cxp 



P— — QO 



2Ai 
A 



M{r) + ^{^ x[i^{R)u2{R)]){r) 



- ^-Kpf 



1 

(31) 



(32) 



We tabulate the Villain potential in the last line before the start of the simulations. Since this 
action contains fewer variables than the previous actions, it is more efficient to run in Monte 
Carlo. 
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4. Numerical Evidence for Gapless Edge 



In order to determine the existence of gapless modes on the edge, we will measure correlators 
of the and 02 variables. 

A technical point: these variables are not translationally invariant because of the way in 
which we introduced fluctuating boundary conditions to enforce a zero total current in the J 
variables. In particular, these variables appear as V^i0 everywhere in the system, except at the 
cuts where we added the 7 variables, where they appear as V^i^ — 7^;. Assuming here that 
takes values 0, 1, . . . , — 1 and the fluctuating boundary condition is between = — 1 
and = 0, we can make a change of variables — >■ 0, 

0(r)=0(r)+^7^r^/L^. (33) 

This gives V^0(r) — 7p(5r^=Lj,-i = Vp0(r) — 7^/Lp, and the action becomes translationally 
invariant along the edge in the variables, so we will use them to measure correlators. 
In the 01 variables, we measure the correlator 

Xi(r- - r') = (e**iWe-''^i(''')). (34) 

Gapless modes exist when xi('' ~ has algebraic decay along the edge. We have taken r' — 
r + m/t, where m is an integer such that Q < m < L and /i is a unit vector in either the y ox z 
directions. When we present numerical data we have averaged over all directions and all sites on 
the edge. We can choose either edge to measure the correlators at; in this work we have measured 
at X = Xii- Since the edge is defined with respect to the dual lattice, it is not obvious where 
on the direct lattice to perform measurements of the edge states. If the sites of the direct lattice 
are indexed by (x, y, z), in our numerics we defined the sites of the dual lattice to be located at 
(a;+l/2, y+1/2, z + 1/2). Using this definition and Eq. ( |25] ) we can see that the edge effects will 
be most noticeable at a; = — 1/2 or x = Xb. — and so we measured the 0i correlators 
ata; = Xfl - 1/2. 

We can gain some insight into the behavior of these correlators by comparing them to spin- 
wave theory with the action 

•S'edge, spin-wave [01 ] = j dydz^ [(Vy0i)^ + (V^^i)^] . 

This is a good approximation to Eq. (IZST i in the limit of small Ai and A2, where vortices in the 
01 variables (which are equivalent to Qi variables) have a very large energy cost to cross the 
edge. Indeed, in the limit A2 ^ one can see that no vortices Qi are crossing the plane where 
we measure the correlators, and the spin wave theory is appropriate on this plane. The spin-wave 
theory predicts the algebraic decay exponent of xi to be 

b^, « l/(2^Ai). (35) 

We expect our data to be consistent with this prediction at small Ai and A2. At larger A2 there 
will be more vortices Qi in our system which increase the decay exponent of the 0i correlations, 
while at larger Ai in the bulk the correlations can also develop through the bulk terms. In our 
numerics taken at Ai ~ A2, the first effect dominates and the extracted power law exponent is 
larger than this spin wave prediction. 
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To detect gapless modes in the J2 variables, we measure the correlator: 



Xi{R - R') = (e*'='^2(K)g-ic02(/?')). (36) 

The c in Eq. ( l36b is the same as the one defining the quantum Hall state with a^y = 2c/ d. As 
discussed after Eq. ( |28] |. in the limit where A2 is very small we can change a 02 variable by 27r/c 
while also redefining p2 and jTi variables and incur only a very small energy cost. Therefore 
for c > 1 if we only measured correlators of e**^^ '■^^ we would not see any order since each (j)2 
variable will be randomly distributed in one of c orientations, and it is the fluctuations around 
these orientations that will be power-law correlated. Thus, our direct analysis suggests that for 
c > 1 only c-tupled "molecular" states of type-2 bosons can propagate along the edge. 

The X2 data can be compared to the following "spin-wave" reasoning. Consider for simplicity 
edges of the a^y ~ 2 quantum Hall state and reformulation Eq. (l3Tt . In this case we can easily 
see the mechanism which forbids vortices on the edge. In the limit of small Ai, we expect the 
currents of M to be zero everywhere away from the edges, while they satisfy a conservation law 
on the edge, J2i_i=y z ^/-t^'^t^ ^ 0- The action on the edge at X = Xji has the structure 

5edge[</'2, Mj,, M,] = ^^^^ [^/^-^^ - 2^P^]2, 

/?€edge;/j.— y,j3 

with Py = and = -My. For such a 2D XY model, the curl of the field iPy,Pz) is the 
vorticity, but this curl is precisely X]^=y z ^p^^/j' which is zero. Thus, in this limit, the action 
is like an XY model with completely prohibited vortices. Away from this limit, we expect the 
vortices to still be effectively prohibited, and the spin wave treatment is now justified. From 
examining Eqs. (ISTl i and ( |28] |, we see that at small Ai and A2 such spin-wave theory predicts the 
exponent 

6^, sa 2T:Xi(f (37) 

for the algebraic decay of X2 (R — R')- 

We now present results for xi and X2- In all the results in this section, we show data on 
the line in the parameter space where Ai = A2 = A. All data was taken with a system size of 
L = 20. We know that we are in the quantum HaU phase when A is small, because it is here 
that the Q variables are gapped. In Sec. |5] we will present the phase diagrams of these models 
and will see that, for all rj ^ 0, the system is in this quantum Hall state for X<fi < 0.33 (the 
precise value depends on c and d, but not very sensitively), while for 77 = the system will be in 
a trivial insulator state in this parameter regime. Therefore at small A the edge we are studying 
is between a trivial insulator and a quantum Hall insulator We have measured the correlation 
functions for Ad'^ = 0.07, 0.1, 0.15, 0.2, 0.25 and seen algebraic decay. [We quote Xd^ because 
at larger d it takes smaller values of A for the Q variables to become gapped, as can be seen in 
Eq. (fTSt.l We plot x on a log-log plot, so that if it decays algebraically we will see straight lines. 
Due to the finite size of the system the lines will not be perfectly straight. We compensate for 
this by replacing R~ R' with the 'chord distance' Djm''- 

L 

Drr' = — sm 

TT 

which is often used in studies of (1h-1)D systems. With this substitution the plots should exhibit 
straight lines, and we found that this works very well in the present cases. 
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Tr\R~R'\ 



(38) 



In Fig.|2]we show xi and X2 for the case where 77 = 1. The plot contains straight lines, which 
means that we do have algebraic decay. We can see that the slope of these lines, and therefore 
the decay exponents, depend on A. In taking the data for this plot we used the reformulation 
given by Eq. dSTI ) for the X2 measurements. This reformulation is more efficient than the one in 
Eq. ( |28] |, and allows us to obtain better statistics. However, we have compared the results of the 
two reformulations and found them to be consistent. Note that as A is increased, xi decays more 
slowly while X2 decays more quickly. This behavior suggests that (pi and (f)2 behave as conjugate 
variables, which is one of the predictions of the A'-matrix theory discussed in [Appendix B| 

Figure [3] shows the same measurements for 77 = 2. Again we see evidence of algebraic decay 
with exponents which depend on A. As above, we used the reformulation in Eq. (l3ll for the 
X2 measurements. Note that for this data it is important to measure correlators of exp(i2(/)2), as 
defined in Eq. ( |36] |. As expected, we found that single-boson correlators decay exponentially in 
this case, and only pair-boson correlators show algebraic decay. 

From the above data we can extract the exponents of the algebraic decay. We fit the above 
data to the function 



with A and b^^ parameters of the fit. We analyzed X2 similarly. Figure |4] shows plots of these 
exponents for 77 = 1 and 77 = 2. The decay exponents for both xi and X2 are slightly above the 
spin-wave prediction at small A (e.g., within 10% for A = 0.07). At large A the fitted exponents 
differ significantly from the naive spin-wave predictions, though within an order of magnitude. 
In [Appendix B| we discuss a phenomenological understanding of the edge which predicts that 
the product of these exponents in the integer quantum Hall case should be equal to 1. We can see 
from Fig.|4|that the products we measured are approximately equal to one. 

Figure [5] shows X2 for 77 = 1/3. The existence of straight lines in this plot implies that 
we have gapless modes in a fractional quantum Hall system. We acquired this data using the 
reformulation in Eq. (|28^ . We were unable to measure xi for the fractional cases because the 
decay exponents were too large. Recall that in order to have gapped Q variables we need XdP < 
0.33; for = 9 this leads to a small A, and the spin-wave theory estimate Eq. (i35T l tells us that 
this leads to large exponents for xi- The presence of the factor d? in the product b^^^h^^ « d?, 
which we obtained here by the direct analysis, is an indirect manifestation of the fractionalization 
when parameter d > 1. Indeed, in [Appendix B[ we use a phenomenological model of the edge 
and the fractionalization of the particles to show that the product of these exponents should 
be equal to d^, while if there is no fractionalization the product will be equal to 1. In our 
lattice model of the edge and the specific parameterizations of the potentials, we have found 
that the values of h^^ are numerically similar to those in the integer case, but the values are 
much larger This implies that the product of the exponents is greater than 1, providing indirect 
evidence for fractionalization. 



5. Phase diagrams 

In previous 

works s mm 

we have studied actions similar to that in Eq. (|2]i. However, in 
those works we considered only the case where 6{k) was equal to a rational constant multiplied 
by 2tt. These are precisely the actions of the Q variables which appeared in Eq. (jSj, and from 
now on we will refer to them as "statistical" actions. The statistical variables are quasiparticles 
in the quantum Hall phases. In Refs. l33Li37i1 we did not attempt to connect the statistical actions 
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Figure 2: Correlation functions a) xi and b) X2, plotted against the chord distance of Eq. | |38K on a log-log scale, for 
rj = 1 edge. Error bars come from comparing runs with different initial conditions. The straight lines imply that we have 
algebraic decay in the correlation functions and therefore the edge is gapless. The slope of these lines varies with A, and 
the extracted exponents are shown in Fig.|4] 
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r] = 2 edge 




chord distance 

Figure 3: Same as Fig. |2] but for rj = 2 edge. Once again, we see evidence for gapless modes. Note tliat in b) we 
measure pair-boson e*^'^^ correlators on the edge, while single-boson e"^^ correlators decay exponentially. 
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?7 = 3 edge 




chord distance 



Figure 5: X2 congelation function for rj = 1/3 edge, wliere c = 1, d = 3. Tlie straight lines indicate algebraic decay 
for this fractional case. We were unable to obtain data for xi because it decayed too quickly. 
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to a physical system, and instead focused on their phase diagrams. Therefore we did not specify 
the "vacuum" of physical variables as it does not affect the dynamics of the phase transitions, but 
it is this vacuum which carries the quantized a^y as can be seen from Eq. ([8]l. Furthermore, we 
can now also specify physical charges of the quasiparticles. 

If we start with a statistical action with rational 9g, we can invert the change of variables 
procedure in Sec. ^ with the correct choice of a, b, c, d to get an action with 9{k) ^ fc^, as 
in Eq. ([T]i, which we will from now on call the "physical" action. This is how we obtained the 
specific potentials in Eqs. (fT4l i and ( fTSl l. We know from the previous works the phase diagrams 
of the statistical actions in terms of the variables Q. We can use the change of variables in this 
work to describe these phase diagrams in terms of the physical variables J'. 

Recall that the specific potentials can be classified by the coefficients c and d. When changing 
from the physical variables J' to the statistical variables Q, we need these coefficients as well as 
a and b, but these coefficients are not independent since they must satisfy ad — be = 1. In 
particular, if we have one solution qq, bo then this constraint tells us that 

a = flo + mc, 

b = bo + md, (40) 

are also solutions if m is an integer The statistical actions can be classified by their statistical 
angle, which for our specific potential choices is given by dg = 2'iTb/d. Therefore each physical 
action can be related to multiple statistical actions by our change of variables. However, these 
statistical actions differ only in that 9g can be different by an integer multiple of 2tt. Such a shift 
will have no effect on the partition sum in the Q variables. Therefore all of the statistical actions 
which can be related to a given physical action have the same behavior. 

We can also see that multiple physical actions whose ratio c/d differ by an integer can also be 
mapped to the same statistical action: that is, the actions in terms of G particles are essentially the 
same except for the "background" quantum Hall conductivity changing by an even integer. 
These physical actions are related to each other by adding an "integer quantum Hall layer" to the 
system without changing the properties of the fractionalized excitations. 

To begin the discussion of the broader phase diagrams of our models, it is useful to consider 
also the action in terms of the dual variables Qi, Q2: 

^ - ^El7^hl2iWI' + ^2|Q2(fc)p]+*E^2i(-^)-«22(fc), (41) 

where Q2 = V x ag2. This action comes from dualizing the J\, Ji variables with the specific 
potentials in Eqs. (fT4] i and (flST l. [Eq. (15) in Ref. fsl!] contains this action with general poten- 
tials.] This action can also be obtained by applying the modular transformation (0, —1, 1, 0) to 
the original action. Note that the d in Eq. (HTI) corresponds to the parameter in Eqs. (UTt -dlSt. 
and is not related to the modular transformation used to obtain this phase. The Q variables are 
vortices with the usual long-range intra-species interactions, vg (fc) ~ in momentum space, 
and we see that these vortices are gapped for large Ai and A2. The statistical angle in the third 
term of Eq. (HTt is a rational number However, this rational number is unique to the specific 
choices we made in Eqs. (fl4l i-(fT5Tl, and small short-range modifications of the model can lead to 
a different statistical angle for the Q variables. This is unlike the rational Qg in the quantum Hall 
insulators which is robust to short-range modifications of the potentials. The difference comes 
from the quahtative difference when applying Eq. (|9]l to generic short-ranged v(k) ~ const and 
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9{k) ^ k"^ in the cases d = (duality to only vortices) and d 7^ (more general modular 
transformation). 



5.7. Models with c/d — n, aly = 2n 

We now turn to detailed descriptions of the phase diagrams. First we discuss the case where 
the conductivity in the quantum Hall phase is quantized as an even integer. In this case we start 
with a physical action with the potentials Eqs. (fT4]i-(fT5]l with parameters c = n, c? = 1 for n 
an integer We can get a statistical action by applying the modular transformation (1, 0, n, 1), 
and this gives a statistical action with dg ^ and a background Hall conductivity of a^y = 2n. 
Therefore in the statistical action we have a system of two uncoupled loops, which is a system 



that is well understood. Il39[ 14011 In this system when A; < 0.3325, the variables Qi are gapped, 
and when Xi is greater than this value the Qi are condensed. All phase transitions are second- 
order XY transitions. Figure |6] shows this phase diagram. Since we know the behavior of the Q 
variables everywhere in the phase diagram, we can now deduce the behavior of the physical J' 
variables. In the lower left corner phase we have seen that the J variables are in a quantum Hall 
phase. 

To understand the rest of the phase diagram we must make more precise our earlier definitions 
of what makes a variable "gapped" or "condensed". When a variable is in a phase in which it 
is gapped, the energy cost for having large loops of that variable becomes arbitrarily large and 
only small loops are present. When a variable is condensed the energy cost for forming loops is 
small. A variable is condensed if and only if the variable dual to it is gapped. In some phases a 
variable will be neither condensed nor gapped in the above sense; instead the variable can be part 
of a composite object which is condensed or gapped. This is the situation for the J' variables in 
the quantum Hall phases. We bring this out to indicate that there are more cases than just given 
by binary choice of J'i variable being gapped or condensed. In all cases, the precise meaning 
is provided by finding appropriate transformation that leads to a description in terms of gapped 
particles only. 

With this in mind, we can interpret the rest of the phase diagram. We begin with the phase 
in the upper-right corner, where Ai and A2 are large. We can see that in this phase the potentials 
seen by the Q variables in Eq. (HTl i become arbitrarily large, so both of these species of variable 
are gapped. Therefore both species of J' variable are condensed and this phase is a superfluid. 
The conductivities ct^J. and ct^^ diverge in this phase, while the Hall conductivity cr^J.^ is non- 
universal. 

We now study the off-diagonal phases in Fig. |6] For simplicity we discuss the phase in the 
lower right comer where Ai is large but A2 is small (in fact Ai can become arbitrarily large in 
this phase). The upper left comer is similar with the indices interchanged. From Eq. (fl4l we 
can see that when Ai 00, V2{k) j^. Since A2 is small in this phase, V2{k) can become 
arbitrarily large and the J'2 variables must be gapped. In addition, we can see from Eq. (fTsT l that 
the Qi variables see an arbitrarily large potential and are therefore gapped, so the jTi variables 
are condensed. From the above results, we can conclude that this phase is a trivial insulator in 
the variables and a superfluid in the J'i variables. 

Finally, we note that the specific model, Eqs. (fT4l i-(fT5]l. with c = n, d = 1 does not realize 
the trivial insulator phase with both J'i and J2 gapped. Of course, we can obtain such a phase 
by different modifications of the potentials, e.g., by adding large repulsive pieces to both vi and 
V2, and it would be interesting to study such models in the future. 
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Figure 6: The phase diagram for the model with the potentials of Eqs. U4t - U5t and c = n, d = 1. In the lower left 
phase the Q variables are gapped and we have the integer quantum Hall phase with cr^^ = 2n. In the upper right phase 
the J7 variables ai'e condensed, and we have a superfluid. In the off-diagonal phases, one of the S variables is condensed 
and the other is gapped. 
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5.2. Models with d ^ 1 

In Ref. Issll we studied a statistical action with 9g = 2tt/3. The phase diagram for this model 
is shown in Fig.|2] In the lower left comer the Q variables are gapped and this is the fractional 
quantum Hall phase. Any physical action with c? = 3 and c = 1 + 3m (for m an integer) can 
be related to this statistical action by our modular transformation; here we will discuss the case 
where c = 1. In this case the change of variables needed to get from the J' variables to the Q 
variables is (0,-1,1,3). The fractional quantum Hall phase will have cr;J.^ = 2 • i and excitations 
carrying respective fractional charges of 1/3 and mutual statistics of 27r/3. 

We know from our previous numerical studv llssl that in the middle phase the variables dual 
to the Q variables are gapped. We can also compute the action for the variables dual to the 
Q variables and see that it has the same potential as the action for the J' variables. The two 
actions also have values of d{k) which differ only by an integer multiple of 2tt. Such difference 
will translate to factors e^'^' in the partition sum and therefore will not contribute. Therefore 
if the variables dual to the Q variables are gapped then the J variables should also be gapped. 
Therefore this middle phase is a trivial insulator in the J' variables. This can be confirmed by 
measuring the conductivity numerically in this phase. 

In the upper right corner phase we can see from Eq. (BTi that the Q variables are gapped, 
and therefore the J' variables are condensed and this phase is a superfluid. In our previous 
work lssll we found that transition between the trivial insulator and the superfluid is a pair of 
XY transitions, while we found more complicated behavior at the fractional quantum Hall-trivial 
insulator transition. 

The structure described in the previous paragraph holds for any set of physical variables 
which can be mapped to a statistical action with dg = I-k jm, with m an integer. One exception 
is m = 2, where our specific model with c = 1,6? = 2 has an additional symmetry in the 
Q variables which prevents the existence of the middle phase, cf. Fig. 1 in Ref. 13711 . This is 
discussed in detail in Refs. 134 , 37], while here we note that generic perturbations to our original 
model will break this symmetry and open a sliver of the trivial phase in the phase diagram. 

Finally, for more complicated fractions cjd, our model will have multiple phases in the mid- 
dle of the phase diagram, resembling hierarchy of phases that we found in U{\) x C/(l) loop 
models with marginally long-ranged interactions and modular invariance. ll34l We expect that the 
"middle" phase at the largest A is a trivial insulator, while the other phases are various quantum 
Hall states. For example, in our model Eqs. (fT4l-(fT5t with c = 2, = 5, we found the following 
sequence of phases upon increasing Ai = A2: fractional quantum Hall insulators cr^.^ — 2-2/5 
and = 2 • 1/2, trivial insulator = 0, and superfluid. It would be interesting to explore 
such phase diagrams and phase transitions in more detail in the future. 



6. Discussion 

In this work we have presented physical [/ ( 1 ) x J7 ( 1 ) bosonic models which realize insulating 
phases with a quantized Hall conductivity that can take both integer and fractional values. In the 
fractional case, we also have excitations carrying fractionalized charges and non-trivial mutual 
statistics. We have shown how to study these models in Monte Carlo and found evidence for 
gapless edge modes. We have also presented broader phase diagrams of our models. 

The action in Eq. ([T]i can be derived from a local Hamiltonian, as shown in Appendix C 



When we included an edge in our action by varying ?7(i?) in Eq. ( fTSl ), we do not know pre- 
cisely how that edge is realized in the physical Hamiltonian. It is possible that this method of 
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Figure 7: The phase diagram for the model with c = 1, d = 3. In the lower left phase the Q variables are gapped and 
we have a fractional quantum Hall phase with tr^y = 2 ■ ^ . In the upper right phase the variables are condensed, 
which implies a superfluid. In the middle phase the J7 variables are gapped and we have a trivial insulator. This figure 
is reproduced from Ref. l33h . but the phases have been re-inteipreted in tenns of the physical vaiiables discussed in this 
work. 



including an edge changes the Hamiltonian near the edge in such a way as to create gapless 
modes which are not due to the bulk topological state of the system on one side. For example, 
a local strengthening of boson hopping along the edge could lead to gapless (1+1)D Luttinger 
liquid modes. Irrespective of the microscopic details, the edges that we studied do produce the 
quantum Hall a^y, so at least some of the observed properties are due to the topology of the 
bulk phases. Including edges using different methods, and confirming that the observed gapless 
properties are not artifacts of the method used in this work, is a possible subject of future re- 
search. There is however some evidence that our gapless modes are due to topological effects. 
In the 7y = 1/3 case we expect topological gapless modes to exist on the edge in the bottom left 
corner of the phase diagram, but not in the middle phase. This is precisely the behavior which 
we have observed, with gapless modes disappearing beyond X(P = 0.35. In addition, we have 
observed gaplessness in both (pi and (p2 variables (for the integer case where both signals could 
be detected), which is what we expect if the gapless modes are topological, and in qualitative 



agreement with the phenomenological iiT-matrix theory of the edge in Appendix B 



Our work allows the numerical study of interacting topological insulator phases, and there- 
fore may be able to address many questions about such phases. For example, we could investigate 
the effect of disorder and other perturbations on the gapless edge states. We could also study tran- 
sitions between different topological phases. In Ref. 13311 we have observed unusual behavior at 
the fractional quantum Hall-trivial insulator transition, which could be studied more closely. In 
addition, our model can reahze transitions between different fractional quantum Hall states, as 
well as transitions between integer quantum Hall states and trivial insulators, which are of recent 
interest. ll2lll22ll More generally, it would be interesting to see what other interacting topological 
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phases can allow unbiased numerical studies. Furthermore, since ideas in the present work do 
not rely on Chern-Simons construction specific t o (2+l) D, they may be more readily extended to 
studies of such phases in higher dimensions. Ill ll4l3ll41 1 
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Appendix A. Formal duality procedure 



In the main text we use a duality transform to make a change of variables, for example 
from the physical J variables to the dual Q variables. In this appendix we review this duality 
explicitly for one loop species. 11271 - 13 311 The original degrees of freedom are conserved integer- 
valued currents J(r) residing on links of a simple 3D cubic lattice; V • J(r) = for any 
r. To be precise, we use periodic boundary conditions and also require vanishing total current, 
Jtot = X^r — 0- define duality mapping as an exact rewriting of the partition sum in 
terms of new integer- valued currents Q(i?) residing on links of a dual lattice and also satisfying 
V • Q(i?) = for any R and Qtot = 0. The derivation of this is as follows. As in Sec. [3] 
of the main text, we seek to integrate out the constrained, integer valued J variables. We can 
implement these constraints using Eqs. (|2T1 |-(|23]|. obtaining the following partition sum: 



Z 



exp -S[]] + i 2^i^(r)[V^0(r) - 2^p^( 

(A.l) 

where S[i] is the original action, now in terms of real valued variables. 

Now consider Q variables such that Q = V x p. Clearly there are multiple values of p which 
give the same Q. Two such p, given by p and fp, are related as follows: 



(A.2) 



Here N{r) is an integer-valued field and are integers. We can divide all possible configura- 
tions of p^ (r) into classes, where two configurations are in the same class if they can satisfy the 
above equation. We can separate the above sum over all p^ (r) into a sum over classes (where 
each distinct class is denoted by a fixed member pji (r)), as well as a sum over the members of 
each class, which comsponds to a sum over N{r) and We can then absorb the sums over 
N{r) and A/^ into the definitions of (f){r) and 7^, which changes the limits on their integration 
to (—00, +00). We can then interpret the integration over these variables as producing delta 
function constraints on the j variables. This gives the partition sum: 



(A.3) 
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Note that due to the U(l) symmetry of the action we can fix (/)(r = 0) = 0, N{r = 0) = 0, and 
so there is no deha function at r = 0. 

If we wish to obtain an action entirely in the Q variables, we can now integrate out the j 
fields. In this work, we study actions of the form: 



^I'^] ^lj2''ik)\J{k)\' + tY,J{k)-A{-k), 



(A.4) 



where v{k) is a potential and A is a fixed gauge field coupled to the J' variables (in the main text 
in different contexts it corresponds to the external gauge field A^^^ or the internal gauge field 
02). With this action, the integrations over the j variables are Gaussian, with basic averages with 
respect to the quadratic piece in Eq. (IA.4l i given by: 



OM(fc)v(fc'))o 



'Jk+k'=0 

v{k) 



fk,ij,fj 



\fk\ 



(A.5) 



where f i 



k.fi 



1 — e""" . We then obtain 



^ 1 [2TTQi~k)+B{-k)]-[27TQ{k) + B{k)] 



(A.6) 



where B = W x A. The relation between Eq. (I A. 4b and Eq. ( IA.6b is what we call "duality map" 
in the main text. 

We can also view the duality procedure in the following way, which is useful for the discus- 



sion of the K-matrix approaches in Appendix B.2 Starting with Eq. (IA.3b we can enforce the 
constraints implemented by the delta functions by writing j as the curl of a real valued field: 
j = (V X a)/27r. We can then perform unconstrained integration over all values of a, and the 
partition sum can be written as: 
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(A.7) 



We see that, compared to our original action SlJ], duality corresponds to replacing integer 
valued J' by real-valued (V x a)/27r, and adding a term which couples the a fields to new 
variables Q. 



Appendix B. Relation to other approaches 

In the main text, we obtained the physics of our models directly using exact transformations. 
Here we point out connections to effective field-theoretic approaches that may be more familiar 
to the readers. Our models can provide rigorous framing and testing grounds for such approaches. 

Appendix B.l. Relation to non-linear sigma models with topological terms 

Let us consider our model in terms of the dual vortex variables Q1.2, which has an action 
given by Eq. (HTT i. As argued in Ref. ll42il . a similar vortex loop model arises in an 0{2) x 0{2) 
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theory with a topological term with = 9q. We see that integer quantum Hall states in our model 
(c = n, d = 1) correspond to 9q ~ 2TTn, as proposed in Ref. lll7ll . We also see that fractional 
"c/d" Quantum Hall states constructed in this paper correspond to 9q ~ 2TTc/d. Our models 
thus provide lattice regularizations of such non-linear sigma models with topological terms. Here 
we want to emphasize that such models do not correspond to a single phase; instead, they can 
have rich phase diagrams, as illustrated in the main text. 

Using our action in terms of the Q1.2 variables [Eq. (HTt ] we can obtain an equivalent for- 
mulation of the original physical current model as 

c,r^ ^ ^ 1^ Ai|[Vxa,i](fc)p + A2|[Vxd,2](fc)P 
J^ofgi, a,2, Ji, J2J — x 2^ -jT— (ti.l) 



2V \h 



fe 



i^^i^ X dgi]{~k) ■ dg2{k) +iY^ ■ (3,1 (fc) + J2{-k) ■ dg2{k) 



The above equation is an intermediate step in the exact duality transformation between Qi,2 and 
Ji_2, where agi.2 are auxiliary real-valued gauge fields encoding conserved real-valued currents 



that appear in the duality transformation, see Appendix A We can interpret these gauge fields 



as mediating the interactions of the physical currents J1.2. The action for the gauge fields has a 
mutual Chern-Simons term but also explicit "mass terms"; the latter make the physical current 
interactions short-ranged, as desired. 

Since the gauge fields are real-valued (non-compact), the above representation implies that 



our model action is unitary, i.e., it can arise as a path integral of a quantum Hamiltonian 11431 14411 
(see also direct demonstration in [Appendix C) . In fact, this is true for arbitrary formal parameter 
c/d — > 77. Note, however, that any such model can still have only rational Quantum Hall phases 
of the type described in the main text, with different rational a^y = 2c' /d' in different regimes. 
Indeed, the available transformations to new gapped variables are modular transformations with 
integer entries (a', b' ,c' ,d') and can only produce gapped Quantum Hall phases with rational 
aly. We can still apply this method to analyze microscopic models with any 77, where it is natural 



to try rational approximants c'/d' to 77 and hence natural to expect rich phase diagram. ll34j - l36ll 
while details require case-by-case study. 

Finally, we note close relation to one of the reformulations in the main text. Since we ulti- 
mately want the action in terms of only the conserved currents J7i and J'2, we can perform the 
integration over the gauge fields ctqi and dq2 in any gauge. For example, we can use the gauge 
V • a = and implement it as follows: We replace (V x d)'^ (V x d)'^ + ^(V • (3)^, perform 
unrestricted integration over the d variables, and take the limit of large ^. We can check that as 
long as the currents J' are conserved everywhere, the unrestricted integration over d gives an 
action independent of ^. (A note of caution: the above statement does not hold when the currents 
have sources and sinks - indeed, boson Green's functions are gauge-dependent.) Taking specific 
value ^ = 1 gives Eq. (1211 1 in the main text. In the partition sum, we integrate independently over 
real-valued fields di^2, and we can check directly that this gives the postulated model without 
the recourse to the dual description. We reiterate that the action Eq. (I2TI 1 is not a gauge theory; 
rather, c5i 2 are some local "oscillator" fields mediating short-ranged interactions of the physical 
currents Ji 2- In [Appendix C| we will show how such an action can arise as a path integral of a 
quantum Hamiltonian with only local degrees of freedom and local interactions. 
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Appendix B.2. Relation to K-matrix theories 

Here we show that our modular transformation analysis can be viewed as a derivation of a 
A'-matrix-like theory, albeit with some non-standard structure in general. The transformation 
consists of 1) duality from Ji to Qi; 2) change of variables Qi = d^i — 6C/2, J2 = cT\ — aQ^ 
[inverse of Eq. ©J; and 3) duality from T\ to Q\. We show all steps starting with the physical 
action and give detailed explanations below: 



+ i 



J2 [^1 ■ ^r* +J2-A. 
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^ [^^1 •7-e1^-(&/3 + a^-') 



First, we do not need to specify the microscopic action 5*8.1-. with short-ranged interactions other 
than that it gives the desired "c/d" phase for some parameters; we carry Ss.r. throughout to better 
display all connections. We also keep track of the external vector potentials A'^2- ^^^h step, 
we show explicitly degrees of freedom that form the partition sum. 

1) We treat the duality from the boson current Ji to vortex current Qi as a reformulation of 
the partition sum replacing jTi — > V x /3/(27r), with real-valued gauge field /?, while keeping the 
information about the integer-valuedness of Ji with the help of new integer-valued current Qi 
as shown in Si (see, e.g.. Appendix A 1. 2) Here we change to new independent currents TitG2, 



which is a valid transformation for (a, h, c, d) forming a modular matrix. 3) Finally, we perform 
formal duality from Fi to Qi as in 1): J^i — > V x 7/ (27r) with real-valued gauge field 7, plus 
new integer- valued current Qi, with the result in 5*3. 

As already mentioned, the Ss.r. part of the model is needed to stabilize the phase with gapped 
CJi 2 particles. Once this is achieved, we can view ^3 as a A'-matrix-like theory in terms of gauge 

fields /3, 7, with the matrix K ~ and charge vectors if = (li 0) and = (0, c) for 

coupling to A°T^^ and A'^^ respectively. Because of the mutual Chern-Simons term for the gauge 
fields /3 and 7, we can ignore Ss.r. and use standard A'-matrix formalism|45J to reproduce the 
result Eq. (fT2l i in the main text for the cr^^. We can also reproduce the mutual statistics of the Qi 
and Q2 quasiparticles Eq. (fTTT i and their charges Eq. (ITjt . but note that we must use the specific 
coupling of Q2 to P and include the direct coupling to A'^^ (if a 7^ 0) to obtain correct results. For 
general h and a, these are non-standard features of the theory in Eq. (IB.2t compared to familiar 
AT-matrix theories, but can be accommodated with proper care. Of course, the recovery of all 
results is expected since the A^-matrix formalism is simply gaussian integration over fields (3 and 

25 



7, while the transformations in the main text carry out such integrations implicitly but exactly 
for our model (also including the short-ranged piece Ss.i-.)- We remark that Eq. (IB.2t becomes 
standard if -matrix theory for modular transformations (a, b, c, d) = (0, —1, 1, d) corresponding 
to the simplest integer and fractional Quantum Hall states with a^y = 2/d. We also note that 
there are other ways to arrive at A'-matrix-like theories and the form Eq. (IB.2t is not unique, but 
our solution in the main text is exact and independent of this. 

While we do not learn new information about the bulk properties of our model from this 
if-matrix formulation, the connection is still inspiring. Thus, the form of the A'-matrix suggests 
presence of two counter-propagating chiral edge modes, i.e., one non-chiral gapless mode on the 
boundary. Our direct analysis and Monte Carlo simulations indeed found power law correlations 
in observables on the boundary. We hope that such tractable models can complement the pow- 
erful if -matrix phenomenology and provide detailed testing grounds for more subtle aspects of 
edge theories. 

Appendix B.3. Phenomenological edge theory 

Let us pursue such a if -matrix approach and compare with the properties of the quan- 
tum Hall edges observed in the main text. Consider first the cr^^ = 2 integer quantum Hall 
state. A convenient modular transformation from the physical bosons to gapped quasiparticles 
is (a, b,c,d) — (0, —1, 1, 1). Application of Eq. ( IB. 21 ) gives a standard form of the if -matrix 
theory. 1451 which then suggests the following edge theory: 

5 = y dxdT-^dr(pid.xtp2 + S'int , (B.3) 

where we use Euclidean space-time and orient the edge along the x-axis. Operator e^'''^ creates 
a Qi quasiparticle carrying unit charge relative to A™', while e"^^ creates a Q2 quasiparticle 
carrying unit charge relative to Aj"*. This suggests that e'"'^'^ contributes to the physical boson 
hi, while e"^^ contributes to 62- [Note that inside the quantum Hall region, the microscopic 
Qa variables are different from the physical J'a variables, and hence ipa are distinct from the 
microscopic boson phase variables <f)a used in the Monte Carlo simulations in the main text. 
Our crude intuition is that the other side of the boundary can absorb the difference since the 
vortices are condensed inside the trivial insulator region, and near the boundary we can write 
schematically ba ^ e"^" . Here we do not attempt a microscopic derivation of the edge theory, but 
rather follow the phenomenological if -matrix formalism. 1191 145111 With the above assumptions. 



we see that at the edge the boson phase fields (pi and (p2 behave as conjugate fields, Igfl similar (up 



to numerical factors) to fields (j) and 9 in the familiar single-mode Luttinger liquid theory. 11461 14711 
This is consistent with our observation that the 61 and 62 power law exponents have opposite 
trends (see Figs.|2]|3] and|4]i, which we discuss further below. 

Note that bi and &2 charge conservation prohibits any cosines of the fields (pi and (p2, and 
the edge is robust. |9|l For simplicity, let us consider harmonic interactions of the form ll45ll 

5i„t = J dxdrj- [Uiid^ipi]^ + U2{d,ip2)''] ■ (B.4) 



We can easily integrate out, say, field ip2 and obtain an action for the field ipi only. 



S^, = / dxdT^ 



U2 
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(B.5) 




or a similar action for the field (^2 only. We then deduce the scaling dimensions of the boson 
fields. 

Hence in this model of the edge, the two scaling dimensions satisfy 

A[6i] A[b2] = \ . (B.7) 

As discussed in the main text, the numerical results are approximately consistent with the above 
equation, cf. Fig.lD There is a slight discrepancy which may be due to the presence of additional 
terms in Eq. ( IB.4b or strong finite size effects. 

We reiterate that we do not have a microscopic justification of the above edge theory and the 
specific choices of the interactions. Nevertheless, the above relation between the scaling dimen- 
sions appears to be close to our numerical results for the model studied in the main text. We sus- 
pect that this may be due to a special time -reversal-like symmetry i ^ ^i,<pi ^ —<Pi,(f>2 02 
in the specific model. While the bulk properties and the edges are robust also without such 
symmetry, the details of the interactions will differ and there will be no such exact relation. Nev- 
ertheless, we expect similar general trends just from the fact that (pi and 02 behave as conjugate 
variables. This remark applies also to all discussions below. 

Consider now the cr^^ = 2n states with n > 2. We take n = 2 as an example. In this 

case, any modular transformation that gives gapped quasiparticles is of the form (a, b, c, d) = 
(1 + mn, m, n, 1) with m an integer; this has a ^ and hence our derivation Eq. (IB. 2) gives 
a somewhat non-standard X-matrix theory. Instead of applying the A'-matrix formalism, we 
can obtain a good physical picture of the edge in this case by bringing together two elementary 
aly = 2 quantum Hall states and allowing boson hopping between the two "layers" labeled / 
and //: 



So = I dxdr 
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1 

dxdr— [drVi+dxip2+ + dr(pi~dx^2-] , (B.9) 



5S 



a=l,2 •' 

= - ^ I dxdr t a cqs{V2lp a-) . (B.ll) 

a=l,2 •' 

Here we have omitted intra- or inter-layer interactions other than the boson tunneling between 
the layers that preserves the U{1) x U{1) symmetry. We have also introduced symmetric and 
antisymmetric combinations of the phase fields in the two layers, 

= (^i^)±¥.i^^))/V2. (B.12) 

We expect that one of the cosines in Eq. (IB. lib is strongly relevant and its coupling will flow 
to large values and will pin the corresponding phase field. Note that since tpi_ and ip2- are 
conjugate variables, only one of them can be pinned. 
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Let us assume that gets pinned; the conjugate variable Lp2- then fluctuates wildly. The 
remaining fields fi/2.+ represent one gapless non-chiral mode, and we can now discuss proper- 
ties of the resulting edge. First, we can write the boson fields hi as 

Thus, we expect power law correlations with the scaling dimension 

1 /[/off 

A[6.] = ^^^, (B.13) 

where are some effective couplings, again assuming interactions like those in Eq. ( IB.4b but 
now for and ip2+- On the other hand, the boson fields 62 contain the wildly fluctuating 
phase (y92- in the exponent and hence the &2 correlations are short-ranged. To obtain power law 
correlations, we need to consider an object carrying two 62 charges: 

(b2f - e*^^"e*^^"' = e*^^^+ , (B.14) 

which has scaling dimension 

A[(M^] = y^. (B.15) 

This provides a phenomenological explanation of our findings in the specific edge model in 
the main text, where single-boson fields bi show power law correlations, while only pair-boson 
fields (62)^ show power law (of course, we can also have situation where the two species are 
interchanged). The above scaling dimensions satisfy 

A[6i] A[(&2)'] = \ ■ (B.16) 

We can compare this result to the second panel of Fig. |4] The results are approximately consis- 
tent, suggesting that our picture of the edge physics is correct and generic. 

We can readily generalize the above argument to the integer quantum Hall case with n > 2. 
We can also see what is happening in the bulk, again starting with decoupled cr^J.^ = 2 "layers." 
As discussed in the main text, in each layer we have a condensate of bound states of a 62 charge 
and a vortex in bi. Borrowing language from layered superconductors (although we assume that 
each layer can talk to all other layers), vortices in each layer are "pancake vortices" and are 
connected by "Josephson vortices" running between each pair of layers. In the absence of the 
boson tunneling, the pancake vortices in the layers are uncorrected and there is no line tension 
for the Josephson vortices. When we introduce tunneling, the Josephson vortices acquire line 
tension and the pancake vortices align near the same location on the 2d plane. On a coarse- 
grained scale where the (finite) collection of layers is viewed as a "fat" 2d system, such a stack 
of pancake vortices represents a single vortex in the fat system (indeed, the boson phase winds 
by the same amount in each layer). Since we have a 62 charge bound to the pancake vortex in bi 
in each layer, we have n such 62 charges bound to this single 61 vortex in the fat system. This 
reproduces our picture of the physical origin of the a^y = 2n integer quantum Hall state as a 
condensate of bound states of n charges and a vortex. 

Turning to the fractional quantum Hall cases, we see that in the cr^^ = 2/d case we can 
use modular transformation (0, —1, 1, d), and Eq. ( IB. 21 ) gives a standard form of the A'-matrix 
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theory. We have an extra factor of d in the edge theory, Eq. (IB. 3b . and now e**^! and e"^^ create 
quasiparticles carrying charges 1/d relative to A^''' and A'^*' respectively, so the microscopic 
boson fields are represented as ha ^ e*'*'^" . Performing calculations similar to the above, we 
conclude in this case 

A[6i] A[62] = ^ . (B.17) 

Finally, we can interpret the = 2c/ d edge by bringing together c more elementary cr^j — 
2/d edges. We again assume that the hi tunneling between the layers dominates and flows to 
strong coupling (although this need not always be the case). We conclude that single-boson hi 
correlations are power law, but only "molecular" c-tupled boson (62)'^ correlations are power law, 
and the scaling dimensions are related by 

A[6i] A[(62r] = ^ . (B.18) 

We were unable to check these relations numerically because the scaling dimensions A[6i] were 
too large to measure. 



Appendix C. Hamiltonian formulation 

Throughout the main text, we worked with the Euclidean action formulation of the model. 
The action for the physical currents is local in (2+l)D space-time and is very convenient for 
analysis. However, it is natural to ask whether this action can be realized as a path integral of a 



local Hamiltonian in 2d. 114811 Below we provide an example of such a Hamiltonian. 

We first specify the physical Hilbert space. Our degrees of freedom reside on two inter- 
penetrating square lattices as shown in Fig. IC.8l We place quantum U(l) rotors on sites r of the 
first square lattice. The rotors are described by 27r-periodic phase variables (/)i(r) and conjugate 
integer number variables ni(r), with commutation relations [0i(r), fii(r')] = iS^r'- We place 
another set of U(l) rotors, described by (/)2(R) and n2(R), on sites R of the second (dual) 
square lattice. Finally, we place harmonic oscillators, described by xi and tt^, on centers of 
links of the first square lattice, which are also centers of links of the second square lattice, e.g., 
i =< r, r + X >=< R, R + y > as illustrated in Fig. IC.8l Here are real-valued coordinate 
variables and tti are conjugate momentum variables, [xe,TTi'] = iSu'- Looking ahead, we will 
use a path integral containing both xe and tt^ . We will view the coordinate variables as fields on 
the links of the first lattice, 

j = X or y, while we will view the conjugate momentum variables as fields on the links of the 
second lattice, 

d2j(R) = Ej-fcTT^^^+k • (C.2) 

Here exy = — e^x = 1 is the 2d antisymmetric tensor and < R, R + j > and < r, r + k > are 
crossing links. 14911 Note that in this Appendix we adopt the following notation: Spatial lattice 
sites are labeled with bold face, e.g., r, R. Spatial directions are labeled with Roman letters, e.g., 
j, k\ space-time directions that appear later will be labeled with Greek letters, e.g., ^, v. Oriented 
fields residing on spatial links are viewed as spatial vectors and are labeled with bold face, e.g., 
ai, a2. 
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Figure C.8: Our Hamiltonian Eq. l lC.3t has U(l) rotors residing on sites r of tlie direct square lattice and U(l) rotors 
residing on sites R of the dual lattice. We also have harmonic oscillators residing on crosses £ of the links of the two 
lattices. The first and second U(l) systems are coupled to the oscillator position and momentum variables respectively as 
if the latter were gauge fields. |49|| With appropriate choices of parameters and additional charge-flux couplings, we can 
induce condensations of bound states of charges and vortices leading to the quantum Hall states discussed in the main 
text. 
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Our Hamiltonian is: 



H — Hfii + Hh2 + Hui + Hu2 + + , (C.3) 
Hhi = -^/iiCos[Vj0i(r) - eiaij(r)] , (C.4) 



Hh2 = -^/i2Cos[Vj02(R)-e2a2j(R)] , (C.5) 
R.J- 

^"1 = 2Z!"i["iW+ffi(VAa2)(r)]' , (C.6) 

r 

^«2 = i^M2[n2(R)+.92(VAai)(R)]' , (C.7) 

R 

Here we introduced various parameters such as boson hopping ampHtudes hi and /i2, on-site 
energies ui and 1*2, and oscillator parameters k and m. The hopping terms couple the boson 
phases and the oscillators as if the latter were "gauge fields." The on-site terms couple the boson 
numbers and appropriate fluxes of the "gauge fields": e.g., flux V A ai = V xCtiy — VySi^r is 
associated with a plaquette of the first lattice, or, equivalently a site R of the dual lattice, and is 
coupled with the boson number ^2 (R) on that site. The corresponding parameters 61,62,51, 52 
will be chosen later Here we emphasize that the model is local in the physical variables (i.e., it 
is not a gauge theory), in the same spirit as Kitaev's toric code model. 

We develop imaginary-time path integral by using Trotter decomposition and insertions of 
unity as follows: 

2 

71-2 {H..T) — -~00 

1 01 (r , t) , 712 (R, t) , (r) ^ {|<?!)i (r , r) , ^2 (R, r) , xe (t 
K,^r+^ = E / ^<^2(R,ri) / Dneir.) 

/ J —TT J —00 

ni (r.Tj^ ) — — oo 
2 

I n-i (r , Ti ) , (?!)2 (R, Ti ) , TTf (ti ) ^ {^7ii (r, r 1 ) , 02 (R, t"i ) , tt^ (r 1 ) 

Here we used one set of variables on "integer" time slices r = int x 6t and conjugate variables 
on "half-integer" time slices ti = t + St/2. We also arranged the Trotter decomposition so 
that the pieces of the Hamiltonian act as c-numbers on the kets of the above insertions of unity. 
Throughout, we omit normalization constants. The remaining inputs to complete the path integral 
formulation in the above variables are overlaps such as 

(X£(T + dr)|7r£(Ti))(7rf(Ti)|x£(T)) = e ^ = g 2. . 1. ^ 

(0i(r,T + (5T)|ni(r,Ti))(ni(r,Ti)|0i(r,T)) = e 2 ^ ^iJi.^ .<Pi 
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and similarly for the second rotor variables. In the above, i =< r, r + j >. 

In the action, we have phase variables (/)i(r,T) residing on sites (r,T) of a (2+l)D cu- 
bic lattice and (t)2 (R, ti ) residing on sites (R, rj. ) of a dual cubic lattice (in the main text, 
such space-time points are labeled simply r and R). We also have boson number variables 
ni{v, Ti) and n2(R, r), which we can view as residing on temporal links of the first and sec- 
ond (dual) cubic lattices respectively and write as temporal components of boson three-currents, 
Jir(r, t) = ni{r, ti) and J72t(R, ti) = 7t.2(R, t + St). We introduce spatial current compo- 
nents using an approach familiar in treatments of XY models; namely, we interpret the cosine 
terms in Hhi and Hfi2 as so-called Villain cosines and write, e.g., 

OO ^2 

Jij (r,r)=-oo 

We can now integrate over the phase degrees of freedom and obtain current conservation condi- 
tions, V • jTi = J2fi=x y T ^ = 0' similarly for the three-current J2- Here and below, 
arrows over symbols denote three-vectors such as J\ ~ [Jixi J\y^ Jit), while bold symbols 
refer to spatial parts such as jTi ~ {Jix,Jiy)- 

We still have the oscillator variables, now labeled ai(r,r) and q:2(R, ti) and residing on 
spatial links of the first and second cubic lattices. At this point, we could also integrate over 
these variables and obtain an action in terms of the boson three-currents only. To facilitate the 
integration and show the connection with the loop models in the main text, we will first write the 
on-site terms by introducing auxiliary fields labeled air and residing on the temporal links 
of the first and second cubic lattices respectively, e.g.: 

/OO ^2 
da2r(R,ri)e-^-*"^-t^^^+^^(^^°^)l . 
-OO 

For brevity, we often omit the lattice coordinates on the fields and imply precise geometric re- 
lation between objects on different lattices; e.g., an oriented plaquette on one lattice is also 
associated with a unique oriented bond on the other lattice crossing this plaquette. 
Putting everything together, the final action takes the form 

S[ai,d2,Ji,J2\ = ^ 

+i^[{'^TOti) AOL2+ .giair(V A a.2) + g2a2r(V A ai)] 

+i ^ [eiJl ■ OLi + JlrO-Xr + 62*72 • OL2 + J2tOL2t\ ■ 

Here the wedge operator is V\ f\V2^ k ^jkVijV2k — vixV2y — viyV2x- By rescahng cii = 
a[/ei and similarly for 0.2, and choosing, e.g., ei = 62 = = l/.g2 = ■\/27r(i/c, we obtain 
essentially the same action as in Eq. (I2TI 1 in the main text. The only difference from the main 
text is that there are additional local current interactions containing hi and /12 couplings, and to 
make the actions identical we only need to take hi and /12 large. In particular, the model is in 
the "c/d" Quantum Hall phase for sufficiently small k and sufficiently large m and large ui, 112- 
Thus, we have provided a Hamiltonian realization for our Quantum Hall phases. 

We can also carry out this derivation when the parameters /ii 2, ui.2, ^^^ m, ei 2, .91,2 vary in 
space; in particular, we can study a boundary between Quantum Hall and trivial insulators. Note 
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that there is significant freedom in how to vary the parameters to achieve different phases even 
within the specific model. For example, we can obtain a trivial insulator by taking ei,2 and gi 2 
to be zero while also taking the hopping amplitudes /ii 2 to be small and on-site potentials ui 2 
large. Alternatively, we can take ei 2 to be very large while gi 2 — > and can reach the trivial 
insulator this way even when the bare hopping amplitudes /ii 2 are large (in this case, the boson 
propagation is scrambled by strong phase coupling to the oscillators). The latter route is closer 
to the model we used in the main text when discussing a boundary between Quantum Hall and 
trivial insulators. Such a boundary model in the present Hamiltonian approach will in general 
differ from that in the main text. Indeed, in the Chern-Simons-like piece for the rescaled d[ 2 
variables, 

r + 5iair(V A ^) + g2a2r(V A -^)] , 

6162 62 ei 

the spatially varying e and g couplings can appear non-trivially under spatial derivatives, and 
the action in general cannot be cast in the form of Eq. ([18). However, we can find a pattern 
of couplings that will reproduce our boundary model in the main text: If we take l/ca — ga — 
^/c/(2jrd) for x G [xaL , XaB.] and 1/ea = (7a = otherwise, and take the region [x2l , X2r] to be 
inside the region [xiL,xi]i],we eliminate terms with non-desired derivatives of the couplings and 
can recast the Chern-Simons-like piece for the rescaled d[ 2 variables into the form of Eq. ([18} 
with ri{R) = c/d inside [x2l , X2b] and zero outside. Thus, we have also provided a Hamiltonian 
realization of the boundary model used in the main text. 

While we universally expect gapless boson correlations on the boundary, detailed aspects can 
be different for different realizations. In this paper, we have focused on the crude demonstration 
of the gaplessness for the specific boundary model in the main text. In future work, it would 
be interesting to examine different realizations and systematically explore all aspects of possible 
edge theories. 
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